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Perturbative NLO and NNLO QCD evolutions of parton distributions are studied, in 
particular in the (very) small- a; region, where they are in very good agreement with 
all recent precision measurements of F^iXjQ^). These predictions turn out to be also 
rather insensitive to the specific choice of the factorization scheme (MS or DIS). A 
characteristic feature of perturbative QCD evolutions is a positive curvature of FJ" which 
increases as x decreases. This perturbatively stable prediction provides a sensitive test 
of the range of validity of perturbative QCD. 

The curvature of DIS structure functions like Q^), i.e., its second derivative with 

respect to the photon's virtuality at fixed values of x, plays a decisive role in probing 
the range of validity of perturbative QCD evolutions of parton distributions in the small- a: 
region. This has been observed recently [U El El [4] and it was demonstrated that NLO(MS) 
evolutions imply a positive curvature which increases as x decreases. Such rather unique 
predictions provide a check of the range of validity of perturbative QCD evolutions. However, 
the curvature is a rather subtle mathematical quantity which a priori may sensitively depend 
on the theoretical (non)perturbative assumptions made for calculating it. Our main purpose 
is to study the dependence and stability of the predicted curvature with respect to a different 
choice of the factorization scheme (DIS versus MS) and to the perturbative order of the 
evolutions by extending the common NLO (2-loop) evolution [3] to the next-to-ncxt-to- 
leading 3-loop order (NNLO) g]. 

The valence = w„, dy and sea w = q, g distributions underlying i^2(x, Q^) are 
parametrized at an input scale Qq = 1.5 GeV^ as follows: 

xqvix^Ql) = Ng^x""- {I - xf"- {1 + Cg^^/^ + dg^x + eg^x^-^) (1) 
xw{x,Ql) = iV^a;''™(l-x)^™(l-f c,„^/J-f d^x) (2) 

and without loss of generality the strange sea is taken to be s — s = 0.5 q. Notice that 
we do not consider sea breaking effects {u ^ d, s ^ s) since the data used, and thus 
our analysis, are not sensitive to such corrections. The normalizations N^^ and Nd^^ are 
fixed by Uydx = 2 and d^dx = 1, respectively, and Ng is fixed via xCS + g)dx = 
1. We have performed all Q^-cvolutions in Mellin n- moment space and used the QCD- 
PEGASUS program for the NNLO evolutions. For dcfiniteness we work in the fixed 
fiavor factorization scheme, rather than in the variable (massless quark) scheme since the 
results for F2 and its curvature remain essentially unchanged p? . 

We have somewhat extended the set of DIS data used in [3j in order to determine the 
remaining parameters at larger values of x and of the valence distributions. The following 
data sets have been used: the small- a; and large-cc [7] HI data; the fixed target 
BCDMS data [8j for F| and using > 20 GeV^ and = Q2(i _ 1) > IQ GeV^ 
cuts, and the proton and deuteron NMC data 0] for > 4 GeV^ and W'^ > 10 GcV^. 
This amounts to a total of 740 data points. The required overall normalization factors of 
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Table 1: Parameter values of the NLO and NNLO QCD fits with the parameters of the 
input distributions referring to (1) and (2). 



the data are 0.98 for BCDMS and 1.0 for NMC. The resuhing parameters of the NLO(MS) 
and NNLO (MS) fits are summarized in Table 1. 

The quantitative difference between the NLO (MS) and NLO(DIS) results turns out to 
be rather small [4]. Therefore we do not consider any further the DIS scheme in NNLO. 
The present more detailed NLO(MS) analysis corresponds to x^/dof = 715.3/720 and the 
results are comparable to our previous ones [3]. Our new NLO (DIS) and NNLO (3- loop) 
fits are also very similar, corresponding to x^/dof = 714.2/720 and 712.0/720, respectively. 
It should be emphasized that the pertur- 
batively stable QCD predictions are in per- 
fect agreement with all recent high-statistics 
measurements of the Q^-dependence of 
i^|'(a;, Q^) in the (very) small-x region. 
Therefore additional model assumptions 
concerning further resummations of sub- 
leading small-x logarithms (see, for exam- 
ple, [in]) are not required [TTlfT^ . 

Figure 1 shows our gluon input distri- 
butions in (1) and Table 1 as obtained 
in our three different fits, as well as their 
evolved shapes at = 4.5 GeV^ in par- 
ticular in the small-x region. In order to 
allow for a consistent comparison in the 
MS scheme, our NLO (DIS) results have 
been transformed to the MS factorization 
scheme. Note, however, that the gluon dis- 
tribution in the DIS scheme is very similar 
to the one obtained in NLO (MS) shown in 
Fig. 1 which holds in particular in the small- 
X region. This agreement becomes even bet- 
ter for increasing values of Q^. This agreement is similar for the sea distributions in the 
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Figure 1: The gluon distributions at the input 
scale Ql^l.h GeV^ and at = 4.5 QeV^. 
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small-x region. Only for x > 0.1 the NLO(DIS) sea density becomes sizeably smaller than 
the NLO(MS) one. The NLO results are rather similar but distinctively different from the 
NNLO ones in the very small- a; region at > Qq. In particular the strong increase of the 
gluon distribution xg{x, Q^) as a; ^ at NLO is somewhat tamed by NNLO 3-loop effects. 

Turning now to the curvature of F2 we first present in Fig. 2 our results for F2(x,Q^) 
at X = 10""*, together with a global fit MRSTOl NLO result |l3], as a function of [I] 

q = logio ( 1 + • (3) 

^" V 0.5 GeVV 

This variable has the advantage that most measurements lie along a straight line [5] as 
indicated by the dotted line in Fig. 2. All our three NLO and NNLO fits give almost the 
same results which are also very similar [3] to the global CTEQ6M NLO fit [TJ] . In contrast 
to all other fits shown in Fig. 2, only the MRSTOl parametrization results in a sizeable 
curvature for i^f . More explicitly the curvature can be directly extracted from 

F^{x, Q^) = ao{x) + ai{x)q + a2{x)q^ . (4) 

The curvature a2{x) = F2{x, Q^) is evaluated by fitting this expression to the predic- 
tions for at fixed values of x to a (kinematically) given interval of q. In Figure 3 
we present a2 (x) which results from experimentally selected g- intervals [H [31 [5j : 



0.7 < (? < 1.4 for 
0.7 <q< 1.2 for 

It should be noticed that the average value 
of q decreases with decreasing x due to the 
kinematically more restricted range ac- 
cessible experimentally. (We deliberately do 
not show the results at the smallest avail- 
able X = 5 X 10~^ where the q-interval is 
too small, 0.6 < q < 0.8, for fixing a2{x) 
in ^ uniquely and where moreover present 
measurements are not yet sufficiently ac- 
curate). Apart from the rather large val- 
ues of a2(x) specific [1 g] for the MRSTOl 
fit, our NLO and NNLO results agree well 
with the experimental curvatures as cal- 
culated and presented in [2] using the HI 
data [6] . Our predictions do not sensitively 
depend on the factorization scheme cho- 
sen (MS or DIS) and are, moreover, per- 
turbative stable with the NNLO 3-loop re- 
sults lying typically below the NLO ones, 
i.e. closer to present data [4]. It should 
be emphasized that the perturbative sta- 
ble evolutions always result in a positive 
curvature which increases as x decreases. 



2 X 10-* < a: < 10^2 

5 X 10-^ < a; < 2 X lO""* . (5) 
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Figure 2: Predictions for ^2(2;, Q^) at a; = 
lO""* plotted versus q defined in ([3|). 

Such unique predictions provide a sensitive 
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test of the range of validity of pcrturbative QCD! This feature is supported by the 
data shown in Fig. 3. Future analyses of present precision measurements in this very 
small-x region (typically 10^^ ^ x < 10~^) should provide additional tests of the the- 
oretical predictions concerning the range of validity of pcrturbative QCD evolutions. 

To conclude, pcrturbative NLO and 
NNLO QCD evolutions of parton distribu- 
tions in the (very) small- a; region are fully 
compatible with all recent high-statistics 
measurements of the Q^-dcpcndcncc of 
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are pcrturbativcly stable and, furthermore, 
are rather insensitive to the factorization 
scheme chosen (MS or DIS). Therefore ad- 
ditional model assumptions concerning fur- 
ther resummations of subleading small-x 
logarithms are not required. A characte- 
ristic feature of pcrturbative QCD evolu- 
tions is a positive curvature a2 (x) which in- 
creases as X decreases (cf. Fig. 3). This 
rather unique and pcrturbativcly stable pre- 
diction plays a decisive role in probing the 
range of validity of pcrturbative QCD evolu- 
tions. Although present data are indicative 
for such a behavior, they are statistically 

insignificant for x < 10^^. Future analyses of present precision measurements in the very 
small-x region should provide a sensitive test of the range of validity of pcrturbative QCD 
and further information concerning the detailed shapes of the giuon and sea distributions 
as well. 
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Figure 3: The curvature 02(2;) as defined in 
(HI for the variable g-intervals in ([5]). 
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